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Abstract 

Active  control  of  the  onset  of  stall  instabilities  in  axial  flow  compressors  is  pursued  us¬ 
ing  bifurcation  analysis  of  a  dynamical  model  proposed  by  Moore  and  Greitzer  (1986).  This 
model  consists  of  three  ordinary  differential  equations  with  state  variables  being  the  mass 
flow  rate,  pressure  rise,  and  the  amplitude  of  the  first  harmonic  mode  of  the  asymmetric 
component  of  the  flow.  The  model  is  found  to  exhibit  a  stationary  (pitchfork)  bifurcation 
at  the  inception  of  stall,  resulting  in  hysteresis.  Using  the  throttle  opening  as  a  control, 
analysis  of  the  linearized  system  at  stall  shows  that  the  critical  mode  (zero  eigenvalue)  is 
unaffected  by  linear  feedback.  Hence,  nonlinear  tools  must  be  used  to  achieve  stabilization. 
A  quadratic  feedback  control  law  using  measurement  of  asymmetric  dynamics  is  proposed 
which  stabilizes  the  bifurcation  and  eliminates  the  undesirable  hysteretic  behavior. 
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1.  Introduction 


* 


Recent  years  have  witnessed  an  increasing  interest  in  axial  flow  compressor  dynamics, 
both  in  terms  of  analysis  of  stall  phenomena  and  their  control  [l]-[6],  [11],  [12].  This  interest 
is  due  to  the  desire  for  increased  performance  which  is  potentially  achievable  in  modern 
gas  turbine  jet  engines  by  operation  near  the  maximum  pressure  rise.  This  increased 
performance  is  at  the  price  of  a  significantly  reduced  stability  margin.  Two  basic  types  of 
instability  are  known  to  occur  in  compression  systems.  One,  surge ,  is  a  low  frequency,  large 
amplitude  oscillation  of  the  mean  mass  flow  rate.  The  second,  rotating  stall ,  corresponds 
to  a  traveling  wave  of  gas  around  the  annulus  of  the  compressor.  This  results  in  very 
inefficient  operation  at  constant  mean  mass  flow  rate  and  pressure  rise.  These  descriptions 
are  for  the  fully  developed  (i.e.,  post-transient)  instabilities. 

Greitzer  [1]  employed  a  nondimensional  fourth-order  compression  system  model  and 
introduced  a  nondimensional  parameter,  B ,  which  he  found  to  be  a  determinant  of  the 
nature  of  post-instability  behavior.  A  global  bifurcation  of  periodic  solutions  and  other 
bifurcations  were  found  for  this  model,  and  were  used  to  explain  the  observed  dependence 
of  the  dynamical  behavior  on  the  parameter  B  [3],  [4]. 

Moore  and  Greitzer  [5]  extended  the  previous  model  to  describe  the  surge  and  rotat¬ 
ing  stall  phenomena  in  axial  flow  compression  systems.  This  model  incorporates  nonax- 
isymmetric  dynamics,  whereas  the  model  of  Greitzer  [1]  reflected  only  axisymmetric  flow 
dynamics,  while  employing  a  nonaxisymmetric  (i.e.,  measured)  steady-state  compressor 
characteristic.  Also  in  [5],  the  general  model  was  specialized  to  a  system  of  three  nonlinear 
ordinary  differential  equations.  The  state  variables  of  this  dynamic  model  are  the  mass 
flow  rate,  pressure  rise,  and  the  amplitude  of  the  first  harmonic  mode  of  the  asymmetric 
component  of  the  flow.  For  the  case  of  a  cubic  axisymmetric  compressor  characteristic, 
Moore  and  Greitzer  found  it  convenient  to  use  the  square  of  the  amplitude  of  the  first  har¬ 
monic  mode  rather  than  the  amplitude  itself  as  a  state  variable.  McCaughan  [6]  performed 
a  bifurcation  analysis  of  this  model,  observing  a  stationary  bifurcation  at  stall  inception. 
This  bifurcation  entails  the  local  emergence  of  a  new  equilibrium  point  from  the  nominal 
one.  This  bifurcated  equilibrium  point  is  not  stable,  and  results  in  a  jump  effect  and  thus 
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hysteresis.  We  note  that  McCaughan  [6]  also  observed  a  bifurcation  to  a  large  amplitude 
periodic  solution  for  the  model  of  [5].  (Although  bifurcations  of  equilibria  are  discussed 
briefly  in  Section  2  and  the  Appendix,  we  refer  the  reader  to  texts  such  as  Iooss  and  Joseph 
[10]  for  a  detailed  treatment.) 

Several  techniques  have  been  proposed  for  active  control  of  stall  instabilities  in  axial 
flow  compressors  (e.g.,  [2],  [11],  [12]).  From  an  analytical  point  of  view,  these  methods 
employ  linear  control  to  delay  the  occurrence  of  stall  or  to  achieve  stall  avoidance.  Of 
course  the  physical  mechanisms  for  controller  implementation  differ  among  the  proposed 
active  control  schemes. 

The  present  paper  begins  with  the  recognition  of  the  importance  of  local  bifurcations 
as  determinants  of  the  nature  of  post-instability  behavior  of  axial  flow  compression  systems. 
The  philosophy  of  the  control  component  of  this  work  is  similar  to  that  of  Abed  and  Fu 
[7].  This  entails  determining  feedback  control  laws  which  ensure  the  stationary  bifurcation 
results  in  only  stable  bifurcated  solutions.  Thus,  even  though  the  nominal  equilibrium  is 
not  stabilizable  within  the  framework  of  linear  theory,  it  may  be  possible  to  stabilize  a 
neighborhood  of  the  nominal  solution  for  a  range  of  parameter  values  including  the  stall 
value  of  the  disturbance  parameter.  It  will  be  seen  that  an  additional  outcome  of  the  control 
laws  proposed  here  is  the  elimination  of  a  hysteresis  loop  which  occurs  in  the  open-loop 
system  model.  The  results  of  this  paper  apply  to  models  more  general  than  that  studied 
by  McCaughan  [6],  in  the  sense  that  the  axisymmetric  part  of  the  steady  state  compressor 
characteristic  is  not  required  to  be  a  cubic  function,  but  rather  an  arbitrary  smooth  single¬ 
valued  function  of  mass  flow  rate.  We  do  employ  a  cubic  model  in  demonstrating  the 
results  for  a  particular  compression  system. 

The  paper  proceeds  as  follows.  In  Section  2,  bifurcation  theory  is  applied  to  study  the 
stability  of  axial  flow  compression  systems  in  the  vicinity  of  the  stall  point.  A  pitchfork 
bifurcation  is  observed  in  the  model  at  the  stall  point.  The  dynamical  behavior  of  the 
compression  system  near  this  bifurcation  point  is  found  to  be  strongly  dependent  on  the 
axisymmetric  compressor  characteristic.  This  dependence  is  exhibited  through  a  formula 
showing  that  the  stability  of  bifurcated  solutions  is  influenced  by  the  derivatives  of  the 
axisymmetric  compressor  characteristic  at  the  bifurcation  point.  In  Section  3,  a  throttle 
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opening  control  law  is  given.  This  control  law  circumvents  the  uncontrollability  of  the 
zero  eigenvalue  of  the  linearized  model  at  stall.  A  purely  quadratic  state  feedback  using 
measurement  of  the  asymmetric  component  of  the  flow  is  given,  and  is  shown  to  result 
in  local  stabilization  of  the  bifurcation  leading  to  stall.  Section  4  contains  bifurcation 
diagrams  for  the  controlled  and  uncontrolled  cases.  Concluding  remarks  are  given  in 
Section  5. 

Notation 

6  -  angle  along  circumference 

Css  -  nondimensional  axisymmetric  compressor  characteristic 

A  -  amplitude  of  the  first  harmonic  of  asymmetric  flow 

tine  -  nondimensional  compressor  mass  flow  rate 

A P  -  nondimensioned  plenum  pressure  rise 

F  -  inverse  function  of  nondimensional  throttle  pressure  rise 

Dxy..h  -  partial  derivative  of  function  h  with  respect  to  the  variables  x,  y, . . . 

(•)'  -  denotes  differentiation 

2.  Bifurcation  Analysis  of  Stall  Inception 


A  third-order  lumped  parameter  model,  in  terms  of  nondimensional  variables,  for  an 
axial  flow  compression  system  has  been  introduced  by  Moore  and  Greitzer  [5].  Using  the 
notation  of  [3],  the  model  is 


dA 

dt 


a 


Css(rhc  +  WA  sin  9)  sin  6d9 


—jf  =  -  A P  +  —  /  Css(mc  +  W A  sin  9)d9 

dt  Z7T  J 0 


dAP 

dt 


1 

452 


{mc-F(7,AP)}, 


(2.1a) 

(2.16) 

(2.1c) 


where  W,  a  >  0  are  two  constants,  7  is  the  control  input,  and  is  associated  with  the  throttle 
opening,  and  B  is  the  nondimensional  parameter  of  the  same  name  introduced  by  Greitzer 
[1],  which  is  proportional  to  rotor  speed.  In  the  sequel,  B  >  0  and  F  is  a  strictly  increasing 
function  with  respect  to  each  of  the  variables  7  and  A  P. 
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Taking  the  compressor  characteristic  Css  to  be  a  cubic  function  of  me,  and  using 
the  squared  amplitude  A2  (instead  of  A,  say)  as  a  state  variable,  Moore  and  Greitzer  [5] 
were  able  to  explicitly  evaluate  the  integrals  appearing  in  (2.1),  giving  a  more  convenient 
dynamic  model.  McCaughan  [6]  performed  a  bifurcation  analysis  of  this  model.  One 
conclusion  of  [6]  is  that  a  stationary  bifurcation  from  the  nominal  equilibrium  point  occurs 
as  the  throttle  opening  parameter  7  is  varied.  However,  a  similar  analysis  for  the  more 
general  case  in  which  the  axisymmetric  compressor  characteristic  is  not  necessarily  a  cubic 
function  of  rile  has  not  been  reported.  In  this  section,  we  extend  the  study  of  [6]  to  the 
case  in  which  the  axisymmetric  compressor  characteristic  is  taken  to  be  a  general  smooth 
function  of  mass  flow  rate.  Note  that  in  Section  4  we  also  employ  a  cubic  axisymmetric 
compressor  characteristic,  by  way  of  illustration. 

Suppose  Caa  is  a  smooth  function  of  me,  and  solve  for  the  equilibrium  points  of  (2.1). 
By  Eq.  (2.1a),  it  is  easy  to  see  that  A  =  0  always  results  in  dA/dt  =  0.  However,  there 
may  be  equilibrium  points  of  (2.1)  for  which  A  ^  0.  First  let  us  consider  equilibrium 
points  in  the  case  A  =  0.  Denote  such  air  equilibrium  point,  the  location  of  which  depends 
on  the  parameter  7,  as  a:0 (7)  =  (0,  ^(7),  AP°(7))7  .  The  values  mQc{ 7)  and  AP°( 7) 
should  then  satisfy  the  relationships  mQc  =  F(y,  AP°)  and  A P°  =  Csa(rnQc).  Under  the 
assumption  A  —  0,  one  such  equilibrium  point  x°(j)  is  referred  to  as  the  unstalled  or 
nominal  equilibrium  point  for  axial  flow  compression  model  (2.1).  Note  that  this  is  the 
normal  operating  point  of  the  system,  the  location  of  which  depends  on  the  throttle  control 
parameter  7.  In  the  following,  we  consider  the  stability  of  the  unstalled  equilibrium  x°(j) 
as  it  depends  on  7.  Viewing  7  as  a  bifurcation  parameter,  we  study  possible  bifurcations 
from  a:0 (7)  at  parameter  values  for  which  stability  is  lost.  These  bifurcations  result  in  new 
equilibrium  points  of  the  model.  At  these  new  equilibria,  A  ^  0,  corresponding  to  stall 
inception. 

Let  X  =  (aq  ,  aq,  x%)T  denote  the  state  variation  of  (2.1)  near  the  unstalled  equilibrium 
point  x°(7),  where  X\  —  A,  aq  —  me  —  m°c{ 7)  and  x3  =  A P  —  AP°(j).  The  linearization 
of  (2.1)  at  x°(7)  gives 


dX 

dt 


=L0X, 


(2.2) 
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with 


/ oC-;,!,,, «.(-,))  o  o  \ 

Lo  =  0  C'„(m%)  -1  .  (2.3) 

V  0  4 h  -^DAPF(~i,isr°)) 

From  (2.2)  and  the  Routh-Hurwitz  stability  criterion,  we  have  the  following  stability 
result. 

Lemma  1.  The  equilibrium  point  a;°(7)  is  asymptotically  stable  for  system  (2.1)  if 
C'ss(i71q(j))  <  0,  while  it  is  unstable  if  C'S3(m°c( 7))  >  0. 

The  foregoing  result  is  useful  in  that  it  classifies  the  nominal  equilibrium  as  being 
stable  or  unstable  depending  on  the  value  of  the  parameter  7.  Since  we  are  interested  in  the 
bahavior  of  the  compression  system  (2.1)  for  values  of  7  for  which  the  nominal  equilibrium 
is  unstable,  we  are  led  to  study  the  behavior  of  (2.1)  for  values  of  7  near  the  critical  value 
at  which  stability  is  first  lost.  An  equilibrium  point  depending  on  a  single  parameter  can 
lose  stability  in  one  of  several  ways.  One  possibility  is  the  actual  disappearance  of  the 
equilibrium  in  what  is  known  as  a  saddle  node  bifurcation.  This  entails  the  merging  of 
the  nominal  equilibrium  with  another  equilibrium  as  the  parameter  approaches  a  value 
for  which  the  system  linearization  becomes  singular.  Although  this  is  the  most  typical 
of  the  so-called  statiiionary  bifurcations,  it  is  not  the  type  encountered  in  this  paper.  In 
this  paper  the  nominal  equilibrium  is  known  to  exist  both  prior  to  and  subsequent  to  its 
loss  of  stability,  Two  bifurcation  mechanisms  can  occur  under  this  circumstance.  One 
mechanism  arises  when  a  real  eigenvalue  goes  from  being  negative  to  being  positive  as 
the  parameter  is  changed.  In  this  case,  the  loss  of  stability  coincides  with  an  eigenvalue 
passing  through  0,  i.e.,  crossing  the  imaginary  axis  at  the  origin.  The  second  mechanism 
occurs  when  a  complex  conjugate  pair  of  simple  eigenvalues  crosses  the  imaginary  axis. 
Although  both  scenarios  result  in  loss  of  stability  of  the  nominal  equilibrium,  they  have 
very  different  implications  as  to  system  behavior  once  local  stability  is  lost.  In  the  former 
case,  a  stationary  bifurcation  generally  occurs,  resulting  in  the  appearance  of  at  least  one 
new  equilibrium  point  from  the  nominal  one  at  criticality.  In  the  latter  case,  an  Andronov- 
Hopf  bifurcation  occurs,  giving  rise  to  small-amplitude  periodic  solutions  near  the  nominal 
equilibrium. 


6 


* 


♦ 


The  Jacobian  matrix  of  (2.1)  at  the  nominal  equilibrium  point  is  given  by  Lq  as  in 
(2.3).  For  this  matrix  to  have  a  pair  of  pure  imaginary  eigenvalues,  it  is  not  difficult  to 
see  that  C'9S{fnPc{ 7))  must  be  positive.  (This  is  seen  by  examining  the  trace  of  the  lower 
right  2x2  submatrix  of  (2.3).)  However,  if  this  were  the  case  then  the  matrix  Lq  would 
have  ctC'ss(rh^y( 7))  as  a  positive  real  eigenvalue,  and  the  equilibrium  would  therefore  be 
unstable.  This  means  that  a  stationary  bifurcation  must  occur  before  an  Andronov- Hop f 
bifurcation  for  the  nominal  equilibrium  a:0  of  the  model  (2.1).  In  the  remainder  of  this 
paper,  therefore,  we  focus  on  stationary  bifurcations  from  x°  and  their  stabilization. 

From  (2.3),  the  linearization  of  (2.1)  has  one  zero  eigenvalue  and  two  stable  eigen¬ 
values  when  C'sa(ihQc( 7))  =  0.  This  implies  a  stationary  bifurcation  may  occur  from  the 
equilibrium  point  x°  for  some  value  of  70.  To  analyze  the  bifurcation  behavior  of  the 
model,  we  employ  a  result  from  bifurcation  analysis  recalled  in  Appendix  A.  This  result 
allows  us  to  derive  conditions  for  existence  and  stability  of  a  bifurcation  from  the  nominal 
equilibrium. 

Let  a;0  be  the  equilibrium  point  at  which  C'ss(m°c( 7))  =  0  for  some  7  =  70.  The 
Taylor  series  expansion  of  (2.1)  at  the  point  (x° ,  70)  is  given  by 

=L0X  +  Q0(X,X)  +  Co(X,X,X)  +  (7  -  7°)Zq  X  +  •  •  •  (2.4) 


Here  L0  is  as  in  (2.3)  and 


Qo(X,X) 


C  0(X,X,X) 


L 1 


aC'f3(mac)x  \x2 
\C's's(m°c)(W2xj  +  2xl) 
-^D(AP)ffi(70,AP°)r 


fC%{m°c)(W*xl+4x1x*) 
C's'l(m°c)aW2xix  2  +  ixl) 
-^D{ap)sF(70,AP0)xI 

'<*fa  C'Jsi^c)  0 

0  hC'J3(m°c) 

0  0 


(2.5) 


(2.6) 


(2.7; 


where 

^  =57F(7°,AP0)  and  02  =  £>Ap7F(7°,AP0).  (2.8) 
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Set  /  =  (1,0,0)  and  r  =  lT ,  the  left  and  right  eigenvectors,  respectively,  corresponding 
to  the  zero  eigenvalue  of  Lq.  The  dynamical  behavior  of  (2.1)  with  respect  to  the  variation 
of  7  near  the  unstalled  point  a;0  is  analyzed  next. 

To  check  the  transversality  condition  lL\r  ^  0  (see  Appendix  A),  we  compute 

lL\r  =a<j>\C"s(m°c).  (2.9) 

The  bifurcation  stability  coefficients  are  calculated,  using  Eqs.  (A.3)-(A.4),  as 

Pi  =  lQo(r,r)  =  0  (2.10) 

and 

P2  =  2l{2Qo(r,6)  +  C0(r,r,r)} 

=^{2DapF(7°,  AP°)(C:.(m%))2  +  (2.11) 

The  next  result  follows  readily  from  Lemma  A.l  and  the  discussions  above.  The  term 
“pitchfork  bifurcation”  used  in  the  theorem  statement  refers  to  a  situation  in  which  two 
new  equilibrium  points  emerge  from  a  given  one,  and  both  occur  locally  on  the  same  side 
of  criticality  (i.e.,  either  for  7  >  70  or  for  7  <  70).  It  is  important  to  note  that  the  stability 
characteristics  of  the  bifurcated  equilibria  determine  whether  or  not  hysteresis  and  jump 
phenomena  result  from  the  bifurcation. 

Theorem  1.  Suppose  C"s{rh°c )  ^  0  and  that  F  is  strictly  increasing  in  each  of  its  variables. 
Suppose  also  that  the  stability  coefficient  P2  given  in  (2.11)  is  nonzero.  Then  the  system 
(2.1)  exhibits  a  pitchfork  bifurcation  with  respect  to  small  variations  of  7  at  the  point 
(x°,7°)  where  C'ss(mQc)  =  0.  Moreover,  if  02  <  0  (resp.  p2  >  0)  the  local  bifurcated 
solutions  near  £°  will  be  asymptotically  stable  (resp.  unstable). 

Note  that  the  formula  (2.11)  shows  that  the  stability  of  bifurcated  solutions  is  influ¬ 
enced  by  the  derivatives  of  the  axisymmetric  compressor  characteristic  at  the  bifurcation 
point.  Theorem  1  quantifies  the  local  dependence  of  post-stall  behavior  of  the  compression 
system  on  the  axisymmetric  compressor  characteristic  (i.e.,  on  the  function  C33).  Bifurca¬ 
tion  diagrams  for  an  axial  flow  compression  system  with  a  particular  form  for  Css  is  given 
in  Section  4. 
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3.  Control  of  Stall  Inception 


From  Lemma  1,  the  unstalled  equilibrium  point  becomes  unstable  after  the  parameter 
7  passes  through  the  critical  value  70.  The  point  (a:0, 70)  is  therefore  called  the  stall  point. 
Moreover,  according  to  Theorem  1  the  local  bifurcated  solutions  near  the  stall  point  might 
not  be  stable.  In  this  case  the  bifurcation  is  said  to  be  subcritical.  If  such  a  condition 
occurs,  the  compression  system  may  exhibit  a  jump  from  the  stable  nominal  equilibrium 
when  the  parameter  7  crosses  the  critical  value  70.  This  results  in  a  hysteresis  loop  in  the 
dynamics  of  the  system  with  respect  to  the  parameter  7  near  the  stall  point.  An  example 
of  such  behavior  will  be  given  in  Section  4. 

In  this  section,  we  seek  feedback  control  laws  governing  the  throttle  setting  near  the 
stall  point  which  prevent  the  occurrence  of  this  hysteresis  or  jump  behavior.  The  control 
design  begins  by  observing  that  jump  behavior  occurs  because  there  is  no  nearby  stable 
equilibrium  after  the  nominal  equilibrium  loses  stability.  A  small-energy  throttle  control 
is  applied  to  the  system,  which  is  equivalent  to  replacing  the  throttle  parameter  7  in 
(2.1)  by  7  +  u  where  u  is  a  control  signal  which  is  of  small-amplitude  near  the  nominal 
equilibrium.  By  examining  the  Jacobian  matrix  Lq  of  the  system  as  given  in  Eq.  (2.3),  we 
find  that  throttle  control  cannot  affect  the  system  eigenvalue  that  passes  through  zero  as 
7  varies  through  the  critical  value  70.  (This  is  simply  the  (1,  l)-element  of  To,  due  to  the 
block  triangular  structure  of  Lq.)  Given  these  circumstances,  we  seek  nonlinear  feedbacks 
which  transform  the  subcritical  bifurcation  at  70  into  a  supercritical  bifurcation,  i.e. ,  to 
give  rise  to  stable  equilibrium  points  locally.  These  equilibrium  points  serve  as  alternative 
steady-states,  so  that  a  jump  to  a  distant  equilibrium  will  no  longer  occur. 

Note  that  such  a  control  law,  although  obtained  through  local  analysis,  mitigates  an 
undesirable  global  effect,  namely  hysteresis.  As  in  Section  2,  denote  by  70  the  critical 
value  of  the  throttle  control  parameter  at  which  stability  of  is  first  lost,  the  design  a 
control  law  ensuring  that  the  local  bifurcated  solutions  near  the  stall  point  are  stable. 
An  important  consequence  of  this  is  that  the  controlled  system  will  not  exhibit  jump  or 
hysteresis  phenomena  near  the  nominal  equilibrium. 

Denote  the  stall  point  by  (x°,7°),  and  let  7  :=  70  +  u ,  where  u  is  the  control  input. 
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Then  we  can  rewrite  system  (2.1)  as  the  nonlinear  control  system 


dA 

a  r2n  . 

dt 

=— — -  /  C33(rhc  +  WA  sin  9)  sin  6d0 
kW  J o 

(3.1a) 

dmc 

1  f2n 

dt 

—  —  AP  +  —  /  Css {dnc  +  WA  sin  6)d0 

2?r  Jo 

(3.16) 

dAP 

1  , 

dt 

~4£2  (mc  +“>AP))' 

(3.1c) 

In  this  section  these  equations  and  their  control  are  considered  without  specifying  param¬ 
eter  values  or  the  form  of  the  steady-state  axisymmetric  compressor  and  throttle  charac¬ 
teristics.  In  Section  4  an  example  is  given  where  these  are  specified. 

It  is  observed  from  (3.1a)  that  A  =  0  is  an  invariant  submanifold  of  (3.1)  regardless 
of  the  choice  of  control  input  u,  implying  that  (3.1)  is  uncontrollable.  Moreover,  it  is  not 
difficult  to  check  that  system  (3.1)  possesses  an  uncontrollable  zero  eigenvalue  at  the  stall 
point.  This  means  that  we  cannot  extend  the  stable  region  of  of  the  nominal  equilibrium 
point  for  a  broader  range  of  values  of  the  parameter  7  by  using  a  linear  state  feedback. 

Next,  we  design  a  control  law  ensuring  that  the  local  bifurcated  solutions  near  the  stall 
point  are  stable.  An  important  consequence  of  this  is  that  the  controlled  system  will  not 
exhibit  jump  or  hysteresis  phenomena  near  the  nominal  equilibrium.  To  make  the  selection 
of  a  control  law  tractable  and  systematic,  we  begin  by  restricting  the  control  to  belong  to 
a  parametrized  family  of  smooth  feedback  control  laws.  Formulas  (A. 3)  and  (A. 4)  show 
that,  in  general,  only  terms  up  to  cubic  order  in  the  state  affect  the  values  of  the  stability 
coefficients  and  jj2  at  a  stationary  bifurcation  point.  This  leads  us  to  limit  the  search 
to  control  laws  containing  only  linear  and  quadratic  terms  in  the  state.  Next,  we  notice 
that  linear  terms  in  the  feedback  control  might  affect  eigenvalues  and  eigenvectors  of  the 
system  linearization  at  the  bifurcation  point.  (However,  recall  the  zero  eigenvalue,  being 
uncontrollable,  would  be  unaffected  by  linear  feedback  terms.)  The  analysis  would  become 
cumbersome,  however,  if  eigenvalues  and  eigenvectors  were  not  fixed,  as  can  be  seen  from 
an  examination  of  formulas  (A. 3)  and  (A. 4).  Thus,  we  are  led  to  seek  a  feedback  control 
law  of  the  following  form: 

u  =q\A 2  +  q2A(mc  -  m°c)  +  q3A(AP  -  A P°)  +  cy4(?fic  -  m°cf 
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+  qs(mc  -  m°c)(AP  -  A P°)  +  qe(AP  -  A P0)2.  (3.2) 

where  the  qi  are  constants.  We  seek  values  of  the  control  gains  qi  such  that,  for  system 
(3.1)  with  u  given  by  (3.2),  02  at  the  equilibrium  z°  for  7  =  7°  will  be  rendered  negative. 

Denote  by  0*  and  0%  the  two  stability  coefficients  of  the  controlled  model  (3.1)  at 
r°  for  for  7  =  7°.  Formulas  (A. 3)  and  (A. 4)  are  also  valid  for  the  controlled  system. 
Using  these  formulas,  it  is  straightforward  to  express  0±  and  0 £  in  terms  of  0\  and  02 ,  the 
stability  coefficients  of  the  uncontrolled  version  of  (3.1).  Indeed,  we  find  that  0*  =■  0\  =  0 
and  02  =  02  +  ^■oiqiD1F('y° ,  AP°)C"s(rh(y).  Note  that  0*  =  0  regardless  of  the  choice  of 
control  gains.  By  Lemma  A.l  in  the  Appendix,  this  implies  that  if  02  0=-  0,  then  a  pitchfork 
bifurcation  must  also  occur  in  the  controlled  system.  From  the  expression  above  for  02 , 
it  is  apparent  that  only  the  quadratic  component  of  the  feedback,  q\A2,  contributes  to 
the  determination  of  the  sign  of  Because  of  this,  we  obtain  a  very  simple  quadratic 
feedback  controller  for  the  system  (3.1),  as  summarized  in  the  next  theorem. 

Theorem  2.  Let  Caa{m°c)  ^  0  at  the  point  (x°,7°).  Then  the  stationary  bifurcation  of 
(3.1)  at  (x°,y°)  can  be  rendered  a  supercritical  pitchfork  bifurcation  by  a  purely  quadratic- 
feedback  control  of  the  form  u  =  q\A2 . 

4.  Illustrative  Example 

In  the  foregoing,  we  have  considered  the  analysis  and  control  of  stationary  bifurcation 
and  hysteresis  for  the  third-order  Moore- Greitzer  model  [5].  In  this  section  we  illustrate 
the  results  for  a  particular  compression  system  model  employing,  for  simplicity,  a  cubic 
axisymmetric  compressor  characteristic.  As  noted  in  Section  2,  the  results  of  this  paper 
are  applicable  to  compression  system  models  with  non-cubic  compressor  characteristics  as 
well. 

In  this  section  we  employ  the  numerical  continuation  and  bifurcation  analysis  package 
AUTO  [8]  to  illustrate  the  compression  system  behavior  for  the  chosen  example.  The 
behavior  of  both  the  uncontrolled  and  controlled  systems  is  considered.  The  compression 
system  example  considered  next  involves  the  cubic  axisymmetric  compressor  characteristic 
given  in  (4.1)  below,  as  well  as  the  throttle  characteristic  given  in  (4.2). 


11 


Let  the  axisymmetric  compressor  characteristic  Cas(mc)  be 

Caa{mc)  =1.56  +  1.5(mc  -  1)  -  0.5(mc  -  l)3.  (4.1) 

Let  F  (the  inverse  of  the  throttle  pressure  rise  map)  be  given  by 

P(7,  AP)  =yy/AP.  (4.2) 

Choose  parameter  values  a  =  0.4114,  W  =  1.0  and  B  —  0.5. 

It  is  easy  to  check  that  there  are  two  values  of  7  for  which  the  compression  system 
(2.1)  has  an  equilibrium  point  (A,  me,  A P)  with  C'a3(rhc )  =  0.  Either  using  the  analytical 
results  of  the  preceding  sections,  or  numerically,  we  can  show  that  these  are  pitchfork 
bifurcations.  One  pitchfork  bifurcation  occurs  for  a  large  value  of  pressure  rise  A P,  while 
the  other  occurs  for  small  pressure  rise.  The  former  pitchfork  bifurcation  is  found  to  be 
subcritical  (i.e.,  it  gives  rise  to  unstable  equilibria)  while  the  latter  is  supercritical  (i.e.,  it 
gives  rise  to  stable  equilibria).  These  observations  are  depicted  in  Figure  1.  Solid  lines  in 
this  figure  indicate  stable  equilibria,  while  dotted  lines  indicate  unstable  equilibria.  From 
this  figure  it  is  clear  that  the  pitchfork  bifurcation  of  most  relevance  in  a  practical  setting, 
i.e.,  the  one  occurring  for  a  larger  AP,  is  subcritical.  Thus  hysteresis  is  expected  to  occur 
in  the  compression  system,  and  a  hysteresis  loop  is  discernible  in  Figure  1. 

Each  of  the  graphs  in  Figure  1  shows  the  effects  of  the  pitchfork  bifurcations  occurring 
in  the  system  on  a  system  state  variable.  Figure  1(c)  is  perhaps  the  most  clear.  From  this 
figure,  we  see  that  a  subcritical  pitchfork  bifurcation  occurs  from  the  nominal  equilibrium 
as  7  is  decreased.  A  second,  supercritical  pitchfork  bifurcation  occurs  at  7  =  0.  Note 
that  Figures  l(a,  6,  d )  also  reflect  these  pitchfork  bifurcations,  but  symmetry  in  the  other 
state  variables  results  in  their  being  double- valued  for  the  bifurcated  solution.  Hence  these 
figures  seem  to  indicate  only  one  bifurcated  equilibrium,  rather  than  two. 

Now,  following  Theorem  2  of  Section  3,  we  choose  the  control  input  u  —  1.0  A2,  which 
is  quadratic  in  the  asymmetric  flow  amplitude  A.  In  Figure  2,  bifurcation  diagrams  for  the 
system  with  this  control  law  in  effect  are  given.  From  this  figure,  it  is  clear  that  hysteresis 
loop  of  the  stable  system  equilibria  of  Figure  1  has  been  eliminated.  The  system  will  tend 
to  operate  at  the  stabilized  bifurcated  equilibria  after  stall  occurs. 


5.  Conclusions 
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The  control  of  axial  flow  compressors  at  the  initiation  of  stall  has  been  studied  using 
the  Moore- Greitzer  model.  The  state  variables  of  this  model  are  pressure  rise,  mass  flow 
rate,  and  the  first  harmonic  of  the  asymmetric  component  of  the  flow.  It  was  found  that 
control  based  solely  on  the  system  linearization  could  not  ensure  stable  operation  past  the 
stall  point,  due  to  uncontrollability  of  the  system  eigenvalue  which  causes  instability  as 
the  throttle  parameter  is  varied.  The  bifurcation  behavior  of  the  model  was  studied  and 
used  as  a  basis  for  nonlinear  control  design.  A  simple  bifurcation  analysis  showed  that 
quadratic  feedback  of  the  first  harmonic  of  the  asymmetric  component  of  the  flow  could 
locally  stabilize  the  bifurcation,  and  hence  eliminate  the  global  hysteresis  loop. 
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Appendix  A.  Bifurcation  Stability  Coefficients 


In  this  appendix  we  recall  a  bifurcation-theoretic  result  on  stability  and  bifurcation 
of  one-parameter  families  of  nonlinear  systems. 

Consider  a  one-parameter  family  of  autonomous  systems 

x  =  f(x,p),  {AX) 

where  x  £  IRn ,  p  £  IR  and  the  vector  field  /  is  assumed  to  be  sufficiently  smooth  in  x  and 
p.  Suppose  that  (A.l)  possesses  a  nominal  equilibrium  point  that  depends  smoothly  on  p 
for  all  values  of  p  in  some  interval  of  interest.  (This  means  we  are  considering  bifurcation 
from  known  solutions  [14].)  For  simplicity  of  notation,  suppose  f(Q,p)  =  0  for  all  p. 
Let  the  Jacobian  matrix  Dxf{ 0,0)  possess  a  simple  zero  eigenvalue  with  all  remaining 
eigenvalues  in  the  open  left  half  complex  plane.  The  parameter  p  is  referred  to  as  the 
bifurcation  parameter. 

The  formulas  and  results  given  below  remain  valid  if  the  nominal  equilibrium  is  not 
the  origin  and  the  critical  parameter  value  is  not  p  —  0;  simply  evaluate  all  quantities 
(including  partial  derivatives)  at  the  adjusted  equilibrium  and  parameter  value. 

The  equilibrium  points  of  system  (A.l)  are  the  solutions  of  f(x,p)  =  0,  and  thus 
depend  on  the  value  of  the  parameter  p.  With  the  assumption  that  Dxf{ 0,0)  is  singular, 
system  (A.l)  may  possess  several  equilibrium  paths  x(p)  emanating  from  the  origin  for 
p  near  0.  Such  a  situation  is  known  as  a  stationary  (or  static)  bifurcation ,  and  the  point 
x  =  0,  p  —  0  is  a  bifurcation  point.  The  solutions  x(p)  ^  0  of  (A.l)  emanating  from  the 
origin  for  p  near  0  are  called  bifurcated  solutions. 

Following  the  notation  of  [7],  we  write  the  Taylor  series  expansion  of  (A.l)  at  (0,  0)  as 

x  =  f(x,p) 

—  Lqx  -fi  Qo(x,x)  -fi  Co(x,x,x)  H - 

-fi  p(Lix  +  Qi(x,x)  -| - )  -fi  p2(L2x  4 - )  -I - ,  (.4.2) 

wher  eQk(x,x)  :=  ^h\D /1kxxf{x,  p),  Cf(x,x,x)  :=  ■^Dflkxxxf(x,  p),  etc.,  are  the  quadratic 
terms,  cubic  terms,  etc.,  of  f(x,p).  Here,  the  quadratic  and  high  order  terms  in  the  ex¬ 
pansion  are  written  as  symmetric  forms.  For  instance,  Qi{x,y)  —  Qu{y,x)  for  each  k  >  0 
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and  any  x,y  €  LRn ■  Moreover,  in  (A. 2)  Lk  :=  j^Dllkxf( 0,0)  for  7*  >  1  and  the  J  acobian 
matrix  Lq  is  assumed  to  have  only  a  simple  zero  eigenvalue  with  the  remaining  eigenvalues 
stable. 

Denote  by  l  and  r  the  left  (row)  and  right  (column)  eigenvectors,  respectively,  of  the 
matrix  Lq  corresponding  to  the  simple  zero  eigenvalue.  Set  the  first  component  of  r  to 
1  (possibly  after  a  reordering  of  state  variables)  and  choose  the  left  eigenvector  l  so  that 
Ir  —  1.  Denote 

A  :=  lQo(r,r)  (A. 3) 

and 

A:=2Z{2<?oM)  +  C0(r,r,r)}  (AA) 

where  6  is  the  unique  solution  to 

L06  =  -Q0(r,r),  (A.  5) 

IS  =  0.  (A.6) 

We  are  now  in  a.  position  to  recall  the  following  stability  criterion  for  system  (A.2). 
This  criterion  applies  only  in  the  case  A  =  0,  and  is  useful  for  studying  pitchfork  bifurca¬ 
tion.  Thus  this  result  is  sufficient  for  the  present  paper.  If  it  happens  that  A  7^  0  then 
another  type  of  stationary  bifurcation,  a  transcritical  bifurcation ,  would  be  expected  to 
occur.  For  details,  see  [7]. 

Lemma  A.l.  Let  the  assumptions  above  hold,  and  assume  the  transversality  condition 
lL\r  0.  Then  (A.2)  undergoes  a  stationary  bifurcation  from  x  =  0  at  p  =  0.  If  A  —  0 
and  A  <  0,  the  bifurcated  equilibrium  points  of  (A.2)  for  p  near  0  are  asymptotically 
stable.  If,  on  the  other  hand,  A  =  0  and  A  >  then  the  bifurcated  solutions  of  (A.2) 
for  p  near  0  are  unstable.  In  either  case,  the  bifurcation  is  a  pitchfork  bifurcation,  and  the 
bifurcated  equilibria  locally  occur  only  for  n(lLir)/ A  <  0. 

Figure  Captions 

Figure  1.  Compression  system  open-loop  operating  points 
Figure  2.  Compression  system  closed-loop  operating  points 
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Figure  2  (rf) 
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